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Abstract 



£NJ ' We consider a general class of models with gauge-mediated supersymmetry breaking in 

which the gravitino is the lightest supersymmetric particle. Several qualitatively different 
scenarios arise for the phenomenology of such models, depending on which superpartner(s) 
decay dominantly to the gravitino. At LEP2, neutralino pair production and slepton pair 
production can lead to a variety of promising discovery signals, which we systematically 
study. We investigate the impact of backgrounds for these signals and show how they can 
be reduced, and outline the effects of model parameter variations on the discovery potential. 
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I. Introduction 



Low-energy supersymmetry (SUSY) can provide a natural solution to the hierarchy prob- 
lem associated with the ratio Mz/M-p^^. If nature is indeed super symmetric, it is important 
to understand the mechanism by which SUSY breaking occurs and is transmitted to the parti- 
cles of the standard model and their superpartners. One possibility is that SUSY is broken at a 
scale ~ 10 11 GeV in a sector which communicates with the particles of the Minimal Supersym- 
metric Standard Model (MSSM) only through gravitational interactions. This has historically 
been the most popular approach, and its phenomenological consequences have been and con- 
tinue to be well-studied. In this paper, we will be concerned instead with a different class of 
"gauge-mediated SUSY breaking" (GMSB) models, in which the messengers of supersymmetry 
breaking are the ordinary gauge interactions . 

Because gauge interactions are flavor-blind, GMSB models are highly predictive with re- 
spect to the form of soft SUSY-breaking interactions. In the minimal model of GMSB j2| ; 
the squark, slepton, neutralino, and chargino masses are determined by only a handful of free 
parameters. The MSSM gaugino mass parameters necessarily have a common complex phase, 
which can be rotated away. Squarks and sleptons with the same electroweak quantum numbers 
are automatically degenerate in mass, up to radiative corrections involving Yukawa couplings 
which can be safely neglected for the sfermions of the first two families. Thus GMSB models 
have the pleasant feature that they are automatically free of excessive non-Standard Model 
flavor-changing neutral currents; this also holds in a large class of extensions and variations of 



the minimal model H-|23j. Furthermore, the sparticle mass pattern is highly constrained even 
in extensions of the minimal model which contain many more parameters. This means that 
sparticle spectroscopy may one day provide for critical tests of GMSB. 

However, the most distinctive phenomenological feature of GMSB models may be that, 
unlike in gravity-mediated SUSY breaking models, the gravitino (G) is generally the lightest 
super symmetric particle (LSP). This is because the scale y/~F associated with dynamical SUSY 
breakdown can be as low as 10 TeV. The spin-3/2 gravitino obtains its mass by the super-Higgs 
mechanism, absorbing the spin-1/2 would-be Goldstino which couples to the divergence of the 
supercurrent with strength 1/F. The resulting gravitino mass is 
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where M = (87rGNcwton) -1 ^ 2 = 2.4 x 10 18 GeV. The next-to-lightest super symmetric particle 



(NLSP) can therefore decay into its standard model partner and a gravitino [24, 25, ^6|. (In this 



paper we assume exact imparity conservation, so that there are no competing decays available 
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for the NLSP.) If the scale yF does not exceed a few thousand TeV, the decays can occur 
within a collider detector volume, possibly with a measurable decay length. Furthermore, even 
supersymmetric particles which are not the NLSP can decay into their standard model partners 
and a gravitino, if no competing decays are kinematically allowed. As we will see below, this 
may be an important consideration. The perhaps surprising relevance of a light gravitino for 
collider physics [24j-|37] can be traced to the fact that the interactions of the longitudinal 



components of the gravitino are the same as that of the Goldstino it has absorbed, and are 



proportional to l/m% (or equivalently to 1/F 2 ) in the light gravitino (small F) limit 



In the GMSB models to be considered in this study, the NLSP is always either a neutralino 
or a charged slepton. In the former case, the lightest neutralino (N\) decays into a photon and 
a gravitino with a width 
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where ki 7 = \Nu cos9\y + N\2 sin#vH is the photino component of N\ (using the notation of 
[38] for the neutralino mixing matrices Nj). The probability that an N% with energy E in the 
lab frame will decay before traveling a distance x is then 

P{x) = 1 - e~ x/L , (3) 



where 



-7 1 ( m N x 




L = 99 x 10 ^ liooSv ) 1^11 E *- - 1 1 ™- (4) 



In principle, one can also have Ni — > ZG or hG, but the corresponding decay widths 
suffer a strong kinematic suppression and can easily be shown to be always negligible within 
the context of the present paper. 

In the rest frame of the decaying Ni, the photon is produced isotropically (independent of 
the spin of Ni) with energy equal to /2. The gravitino still escapes the detector, carrying 
away missing energy. Therefore SUSY discovery signals at colliders involve up to two energetic 



photons and missing (transverse) energy in GMSB models with a neutralino NLSP [|2J]-[36]. At 
the Tevatron, the largest production cross sections typically involve chargino (Cj) and neutralino 
{Ni) production, especially pp — > C^C^ and C x iVjj. One can therefore detect supersymmetry 
using an inclusive 77 + + A signal, in addition to channels with lepton(s) + jet(s) + or 1 
photon. The discovery signatures for SUSY with a prompt decay N% — ► 7G are so spectacular 
that it is possible to set quite significant bounds even with existing Tevatron data. For example, 
in [30 1 it was argued that with the present ~ 100 pb _1 of Tevatron data, it should be possible 
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to exclude chargino masses up to about 125 GeV and neutralino masses up to about 70 GeV 
in a large class of models with a light gravitino, including GMSB models, as long as the decay 
N% — > jG occurs within the detector .0 In ref. pHR , a significant reach was found also in the 
and 1 photon channels. 

In discussing exclusion possibilities at the Tevatron, we must mention that a single unusual 
event |}{| of this general type has been observed at CDF. This event has large (> 50 GeV) 
Et and two central (|r/| < 1), energetic {Et > 30 GeV) photons and two energetic leptons. 
Events with these characteristics are reputed to have very small Standard Model and detector 
backgrounds, and it was pointed out in that this event might be explained by GMSB 

models (and other models with a light gravitino) in terms of selectron pair production. However, 
at least in the simplest types of GMSB models, this interpretation is now perhaps somewhat 
disfavored, since one might typically expect many accompanying events in other channels p9| , 
30, [3l], 35], which have not turned up in recent searches by CDF |4(| and by DO [41]. Perhaps 



a more plausible explanation of the event within the GMSB framework is that it was due 
to chargino pair production pp — > Cy^Cf , as proposed in |30|, |3l[| . Each of the charginos 
can decay either hadronically into qq'^/G or leptonically into £ujG. The latter possibility 
can be significantly enhanced if sneutrinos are not too heavy (although it still seems somewhat 
problematic to explain the kinematics of the observed event). GMSB models can be constructed 
with a neutralino NLSP and a large leptonic branching fraction for C\ (with > mo), but 
not with a minimal messenger sector. The signal for CfCi production in this case can then be 
£ + 1'~^-^Et- Note that in this chargino interpretation of the event, the leptons need not have 
the same flavor. In any case, more data at the Tevatron and at LEP2 will help to test these 
speculations. 

In other GMSB models, one finds that the NLSP is a stau. Here, one should distinguish 
between several qualitatively distinct situations. If tan/? (the ratio of Higgs expectation values 
(H®)/(H®)) is not too large, the lightest stau eigenstate f% is predominantly right-handed and 
is nearly degenerate in mass with the other right-handed sleptons. In the case that ~ 
rrijx R < mf x + m T , one finds that the e R and fx R cannot have three-body decays into f\ without 
violating lepton flavor conservation. Since lepton flavor-changing interactions are automatically 
very highly suppressed in GMSB models with i?-parity conservation and decays through an off- 
shell tau are insignificant, each of the right-handed sleptons decays only into the corresponding 



lepton + gravitino [26], and fi, jl R , e R act effectively as co-NLSPs. (An exception occurs if 

mfi < rrift < mg R , as discussed below.) In this case, all supersymmetric decay chains will 

terminate in fx — > tG or e R — * eG or fl R — > (iG. The formulas for the relevant decay widths 

of slepton into lepton + gravitino are given by simply replacing Wft —* my, — * Ey and 

A quite similar bound on mgj in a large class of GMSB models may be obtained using the recent LEP runs 
at \fs = 161, 172 GeV, as we will see in section III. 
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Ki 7 — > 1 in (||) and (||), so that the decay length in the lab frame for a slepton with energy Ej 
is 



L — 9.9 x 10" ^L-. j - 1) cm. (5) 

At LEP2 energies, the primary discovery process often (but certainly not always, as we shall 
see) involves simple pair production of the NLSP. In that case, E^ in eq. @ and Ej in eq. (|5|) 
can simply be replaced by the electron beam energy at LEP2. 

Conversely, if f\ is the NLSP and tan j3 exceeds 4 to 8 (depending on the other model 
parameters), one finds in GMSB models that is small enough that the decays £r — > f\rl are 
kinematically allowed for £ = e,fx. These three-body decays are mediated by virtual neutralinos 
and are typically not dynamically suppressed, because the bino content of N\ is significant. 
However, we have checked that they can be quite strongly suppressed by phase space, so that 
it is possible for £r — » £G to dominate even if mj — — m T — mi is a few GeV, if y/F is 
not too large and £r — > £N\ is kinematically disallowed. Barring these circumstances, f% acts 
as the sole NLSP, and all super symmetric decay chains will terminate in f± — > tG |34 



An important exception to the preceding discussion occurs if \rrijj — rrifi\ < m T and 
< mj R + Tnt for I = e, fx. Then each of the decays Ni — > jG and fi — > tG have no 
significant competition, and Ni and t\ act effectively as co-NLSPs. 

To summarize, there are four qualitatively distinct scenarios for the termination of super- 
symmetric decay chains in GMSB models. By a slight abuse of language, we refer to these as 
"neutralino NLSP" , "stau NLSP" , "slepton co-NLSP" , and "neutralino-stau co-NLSP" scenar- 
ios, according to whether A only, B only, B and C, or A and B of the decays 

A) JVi -► 7G (6) 

B) n - tG (7) 

C) £ R ^£G (£ = e,n) (8) 

do not suffer competition. The four possible scenarios correspond nominally to the mass or- 
derings (in addition to < m~ for £ = e,fM, which turns out to be always satisfied in the 
GMSB parameter space we consider): 

neutralino NLSP 
stau NLSP 
slepton co-NLSP 
neutralino-stau co-NLSP 



rrift < mfj — m T (9) 

mfj < Mm[m^,mj^] — m T (10) 

fnj R < Minfm^ , mf x + m T ] (11) 

\m fl - m^J < m r ; < mj o , (12) 
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where we have neglected the masses of the electron and the muon. We should note that the 
condition eq. (|l0|) for the stau NLSP scenario is necessary but not quite sufficient, since as 
we have already noted, the decay £r — ► IG can dominate over £r — ► fyrl when the latter 
is kinematically open but suppressed. Thus some models which obey eq. (|l0|) may actually 
behave as slepton co-NLSP models, depending on yF. We have checked that two-body decays 
t — > £N\ always dominate over decays into the gravitino as long as y/F > 10 TeV and the mass 
difference mj— — mi is more than of order 10 MeV for £ = e (and much less for £ = fi, r). 
A similar statement holds for two-body decays N\ — > ££. Hence we will consider only the four 
main scenarios listed above. Other "borderline" cases with small mass differences will have 
similar phenomenology to the cases we do treat. 

In section II of this paper we describe the framework for an exploration of the parameter 
space of GMSB models, with some simplifying but hopefully not overly restrictive assumptions. 
In section III, we give some conditions on the parameters for each of the four different NLSP 
scenarios, and study the possible signals and backgrounds which arise at LEP2 in each case. 
We will mostly consider an option with \fs = 190 GeV and 300 pb -1 per detector. Section IV 
contains some concluding remarks. 



II. Models of gauge-mediated SUSY breaking 



In this paper we will consider the following class of GMSB models. The ultimate source 
of SUSY breaking is parameterized by a gauge-singlet chiral superfield S whose scalar and 
auxiliary components are both assumed to acquire vacuum expectation values (VEVs), denoted 
S and F$ respectively. The superfield S couples to a "messenger sector" consisting of chiral 
superfields <£j, &i which transform as a vector-like representation of SU(3)c X SU(2)l x [7(l)y. 
The messenger sector couples to the SUSY-breaking sector through the superpotential 

W = A i 5# i ¥ i . (13) 

This implies that the fermionic messengers acquire a Dirac mass XiS, while their scalar partners 
obtain (mass) 2 = | *S*| 2 ± |AjFs|. The ordinary gauge interactions then transmit this SUSY 
violation to the MSSM fields, with computable superpartner masses 0, @. Contributions to 
gaugino masses due to each messenger pair $j, $j arise at one loop and are given at the scale(s) 
Q = AjS by 

AM a = ^An a (i)g( Xl ) (a = 1, 2, 3) (14) 

where 

A = F S /S, (15) 
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g(x) = i[(l + x) ln(l + x) + (1 - x) ln(l - a?)] , (16) 

and each 

xi = \F s /\iS 2 \. (17) 

The latter quantities must satisfy < Xj < 1 (so that the lightest messenger scalar does not 
acquire a VEV). Here n a {i) is the Dynkin index for the messenger pair <3?j in a normalization 
where n a = 1 for N+N of SU(N). We always use a Grand Unified Theory (GUT) normalization 
for a\ so that n\ = |i^ 2 for each messenger pair with weak hypercharge Y = Qem — T3. In 
the limit of small X{ and when consist of a 5 + 5 of the global SU(5) which contains 

SU(3)c x SU(2)l x U(1)y, one recovers the results for the original minimal GMSB models 
with ^2n a = 1 and small Xj 0, since g(0) = 1. The function g(x) is always slightly greater 
than 1, but never exceeds 1.044 when x < 0.5, and reaches a maximum of 1.386 at x = 1 [||. 

Contributions from each messenger pair to the (mass) 2 terms of the MSSM scalars 

arise at two-loop order and are given at the scale(s) Q = XiS by 

2 



Am 2 = 2A 2 W^ C a n a (i)f( Xi ) (18) 



where [|| 



ln(l + x) - 2Li 2 (x/[l + x]) + ^Li 2 (2x/[1 + x\) 



+ (x^ -x) (19) 



2 

and C a is the quadratic Casimir invariant of the MSSM scalar field in question, in a normal- 
ization where C3 = 4/3 for color triplets, C2 = 3/4 for SU(2)l doublets, and C\ = \Y 2 . For 
small Xj, one has f(xi) ~ 1, so that again the results of the original minimal model with small 
Xi Q are recovered. 

In order to have a manageable parameter space for our study, we now make some simplifying 
assumptions. First, we consider (except when explicitly noted) only models for which the total 
Dynkin indices of the messenger sector for each gauge group are equal and do not exceed 4: 



n 



Y, «i(») = E = E n ^i) = 1, 2, 3, or 4. (20) 



This assumption ensures that the apparent near-unification of perturbative gauge couplings 
near My ~ 2 x 10 16 GeV is maintained. (Possibilities which do not embrace this assumption 
are discussed in M.) We will also take all of the couplings Aj to be equal to a common value 
A, even though no symmetry can enforce this; variations in the individual Aj only affect the 
MSSM sparticle spectrum logarithmically. This in turn forces all of the x, to be equal to a 
single parameter x, which as a practical matter we require to satisfy 

0.01<x<0.9. (21) 
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With these assumptions, MSSM phenomenology is determined by just 6 parameters 



A, n, x, A, tan f3, and sign(/i) . (22) 
The expressions for the sum of contributions to gaugino and sfermion masses now simplify to 



M a = nAg(x)^ (23) 

47T 



m 



~'-2nA 2 /(x)E(^) 2 Ca (24) 



at the single messenger scale Q meS s = A/x. Equations ( f23[ ) and ( p4| ) are taken as boundary 
conditions for renormalization group (RG) evolution of the MSSM parameters. At the same 
scale Qmessi the running trilinear scalar couplings of the MSSM are taken to vanish (they 
actually receive contributions at two- loop order which are negligible in the first approximation) . 
However, we do not assume that B\x is close to zero at the messenger scale, since it seems likely 
that a mechanism for generating fj, can also generate B[i |3|, ra. We then evolve all of the 
couplings of the MSSM from Qmess down to the electroweak scale, where the parameters Bfi 
and are determined by requiring correct electroweak symmetry breaking. Note that in this 
parameterization, the sparticle spectrum (with the exception of the gravitino mass!) does not 
depend on A, and depends on x only logarithmically through Qmcss and g(x). Of course, allowing 
different Xi,\i would be more realistic and would enlarge the parameter space. However, the 
features of the enlarged parameter space obtained in this way do not differ dramatically from 
the one we consider. The effect of finite Xi in (0) and (|l|) is simply to raise the gaugino mass 
parameters by up to about 25% with respect to the sfermion masses, since g{x)/ \J f(x) varies 
between 1 (for 1) and 1.25 (for x = 0.9). Choosing a specific value of x can be thought of 
as simply parameterizing our ignorance of these effects within a simplified framework. In the 
present paper, the chargino, squark, and gluino masses have no direct relevance, so that the 
practical effect of increasing x is essentially just to lower the slepton masses compared to the 
mass of the lightest neutralino. 

The Goldstino decay constant in this parameterization is given by 

p _^_A2 
XjXj x\ 



This way of expressing F$ is useful because A is relatively well-known from (|23|) and ([24|), 
since it is correlated strongly with the mass of the NLSP which we must presume lies in the 
50 to 100 GeV range in order for SUSY to be relevant at LEP2. An easy estimate then 
shows that the relevant range of A for this paper is from about 10 TeV to about 100 TeV. 
While the dimensionless couplings x and A can be arbitrarily small, they can be bounded from 



S 



above. Furthermore, F$ can be smaller than the full SUSY-breaking order parameter F of 
the complete theory. (In models where F$ arises directly from an O'Raifeartaigh mechanism 
one expects F$ = F, but in models where F$ itself arises radiatively as in one may find 
Fs <C F.) Therefore eq. ( p5[ ) puts a lower limit on F, or equivalently a lower limit on tHq, 
corresponding to a lower limit on the decay length of the NLSP according to eqs. (Q) and (|5|). 
In particular, we note that in viable GMSB models one must have 



> 2 x 1(T 2 eV, (26) 



G ~ 

based only on Aj <, 1 and ttt-nlsp > 50 GeV; this is consistent with cosmological bounds 
on rriQ. For any given sparticle spectrum specified by the parameters A, n, x, tan /?, sign(/i), one 
obtains a lower limit on the NLSP decay length by assuming A <, 1. By taking smaller A with 
A and x held constant one can essentially arbitrarily increase the NLSP decay length while 
holding all other features of the MSSM sparticle spectrum fixed. The NLSP decay length will 
also be increased if Fs < F, as long as a pseudo-Goldstino field which is not absorbed by the 
gravitino (and which is predominantly the fermionic component of S) acquires a large mass. 
In these cases the NLSP decay length L can be made so large that NLSP decays always occur 
outside the detector. 

The statement that A, or more generally the distinct couplings Aj, should be bounded 
from above can be motivated as follows. An estimate of the maximum value of the couplings 
Aj should roughly correspond to an infrared quasi-fixed point of the RG equations. Let us 
consider, for example, a "minimal" messenger sector in which <I>j, $j consist of a 5 + 5 of the 
global SU(5) group which contains the MSSM gauge group. In that case one has couplings A2 
and A3 of S to the 5C/(2)i-doublet and 5f7(3)c-triplet messenger fields, respectively. These 
couplings satisfy the one-loop RG equations 



9 d\2 

1671 -JT = A 2 
at 

16tt 2 — i = A3 
at 



4A|. + 3A| + . . . - 3g 2 2 - ^g\ 
2A 2 . + 5A! + . ..-^gl-± g j 



(27) 
(28) 



Here the ellipses represent the effects of other dimensionless couplings involving S, <I>j, or <3?j. 
These are of course highly model-dependent but will contribute positively to the one-loop (3 
functions, thus only reducing the quasi-fixed point values of the couplings. One can then 
estimate the maximum quasi-fixed point values for (A2, A3) by taking A3 = A2 = \u to be large 
at the putative unification scale and evolving down to the messenger scale. For example if the 
messenger scale is Qmess = 100 TeV, then taking \jj > 2 yields (A2,As) ~ (0.7, 1.1) at Qmess- 
Even if we abandon GUT boundary conditions in this example and choose A3 <C A2 or A2 -C A3 
at the "unification" scale, the maximum values at the messenger scale of the dominant coupling 
are found to be A2 ~ 1.0 or A3 « 1.2, respectively. Of course, S need not be a fundamental 
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chiral superfield, but the fixed point values for Aj should roughly correspond to the maximum 
values, at least in a perturbative effective field theory description. Note that models with 
more messenger fields should respect the rough result A.; < 1. Adding more couplings to 
the mix will effectively only give positive contributions to the corresponding (3 functions when 
compared with the example given above, leading to smaller values for the Aj at the messenger 
scale. A larger messenger sector will cause the gauge couplings to be larger above the messenger 
scale, indirectly resulting in a decrease in the (3 functions for the Aj. However, this can only 
slightly increase the quasi-fixed point behavior of the largest coupling(s), which is determined 
predominantly by what happens near the messenger scale. Therefore we expect that the lower 
bound on the NLSP decay length at LEP to be inferred from A ^ 1 (for a given x) should be 
robust. In particular, we see from (p5|) that for models with a given A, the minimum possible 
decay length L for the NLSP is given by replacing y/~F —> A in (Q) or (|5|). 

In order to understand the parameter space of the GMSB models we have chosen to study, 
we have used a computer program to generate several tens of thousands of models for each of 



n = 1, 2, 3, 4 and random values for the other free parameters in eq. (22). The program proceeds 
by an iterative method that sets the weak-scale gauge couplings and masses, evolves the RG 
equations to the messenger scale, sets the messenger scale boundary conditions, evolves the RG 
equations with associated decouplings at each sparticle threshold back to the weak scale, then 
iterates to convergence (about 4 iterations are typically necessary). Two-loop RG equations 
are used for the gauge couplings, third generation Yukawa couplings, and gaugino soft masses, 
while one- loop RG equations are used for the other soft masses and scalar trilinear couplings. 
Electroweak symmetry breaking is enforced at the scale Q = yjm^ fni 2 , allowing the evaluation 
of \fj,\ and B[i from the Higgs soft masses, tan/3, Mz, and one-loop corrections. This evaluation 
utilizes the one-loop effective potential, which includes the corrections from stops, sbottoms, 
and staus consistently with the one-loop evaluation of the Higgs masses. 

Because this paper is devoted to possible discovery signals at LEP2, we consider only 
models with NLSP mass less than 100 GeV. The lightest supersymmetric particle is always a 
neutralino or a stau throughout this parameter space.Q However, as explained in the Introduc- 
tion, more than one superpartner can act effectively as the NLSP. We find significant regions 
of parameter space in which each of the four scenarios is indeed realized. In the next section 
we will describe in turn some relevant features of the parameter space, including conditions on 
the parameters n, A, and tan /?, for the four NLSP scenarios. In each case we study how the 
SUSY discovery signals may manifest themselves at LEP2. 

2 We find one exception: it is possible to construct models with n = 3 or 4 which have a tau sneutrino NLSP, 
but only if mp < 54 GeV and m~ < 57 GeV and m~ < 95 GeV. We neglect this possibility in the following, 
although a complete analysis which might exclude these models has not yet been performed to our knowledge. 
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III. Signals at the CERN LEP2 collider 



A. The neutralino NLSP scenario 



If equation (g) is satisfied, then the lightest neutralino is the NLSP, and essentially all 
supersymmetric decay chains will terminate in N± — ► jG. If this decay always occurs outside of 
the detector, then collider signatures and search strategies are the same as in the well-studied 



neutralino LSP scenario (see for example [44]). In that case, the only impact of the gauge- 



mediation mechanism for collider phenomenology is that the pattern of supersymmetric masses 
and other soft supersymmetry breaking parameters is restricted in significant ways. 

If the decay N\ — ► 7G occurs within the detector an appreciable fraction of the time, 
then AqAq production can lead to a discovery signal at LEP2. A crucial quantity is then the 
decay length L in equation (Q) with E = yfs/2. For L greater than a few centimeters, the 
LEP detectors may be able to resolve the distance from the interaction point to the N\ decay 
vertex where the photon originates. If this can be done reliably, there should be essentially 
no backgrounds to the signal. However, as L increases, a larger fraction of events will occur 
outside the detector, decreasing the efficiency accordingly. Of course, any analysis of this 
situation would be highly detector-dependent, as the different LEP detectors have varying 
geometries and photon direction resolution capabilities. For a rough study we suppose that 
photons resulting from decays N\ — ► 7G can be detected if they occur within 1 meter of the 
interaction point. Using eqs. (||),(||), we can then estimate, as a function of the parameter L, 
the probabilities that both, or only one, of the photons can be detected in each event. These 
probabilities are shown in Fig. [I] for 10 -3 cm < L < 10 5 cm. (This range corresponds roughly 
to 1 > xX > 10~ 4 for the models described in section II with < 90 GeV and y/s = 190 
GeV.) Note that the probability to observe one of the two photons in each event may exceed 
0.1 even for L greater than 10 meters. Remarkably, an observable signal with displaced (not 
originating from the interaction point) single photons might occur for L up to several tens of 
meters, depending on the AqAq production cross section (typically in the tens or hundreds 
of fb), the integrated luminosity achieved, and the specific detector being used. For L in the 
several centimeter to several meter range, one can hope to observe both one photon and two 
photon events with displaced vertices. 

For the remainder of this subsection, we assume that L is less than a few tens of centimeters, 
so that essentially all decays Aq — > 7G occur within the detector. We consider only models 
for which < 100 GeV (so that AqAq production can be possible at LEP2) and > 70 
GeV (motivated |30| by the non-observation of excessive 77 + fx + V events at the Tevatron 
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Figure 1: The probability that exactly one (two) of the photons from the decay Ni — ► 7G 
originates within 1 meter of the AqAq production vertex is shown as the dashed (solid) curve, 
as a function of the decay length L. 



[40, [Hfl). Then we find that within the framework of GMSB models assumed here, neutralino 
NLSP models can only be constructed if 



n 
n 
n 



40TeV < A< 80TeV 
24TeV < A< 40 TeV 
20TeV < A< 24 TeV 



tan/3 < 35, or 
tan/3 < 18, or 
tan/3 < 10 and \i > 0. 



(29) 
(30) 
(31) 



These conditions are necessary, but not sufficient. The upper and lower bounds on A correspond 
to those onniju in the obvious way. The upper limits on tan (3 follow from the requirement that 



+ m T ; for larger values of tan/3 the mixing in the stau (mass) 2 matrix becomes too 
large to allow this. In all cases, one finds that \fi\ > M%, so that the NLSP has a significant 
photino component, with 0.4 < < 0.85 in all models of this type, and ki 7 > 0.55 in 
the minimal model (meaning n = 1). We have verified that for all neutralino NLSP models 
accessible at LEP2, the BR(A r i — > 7G) is in practice indistinguishable from 100%. In the models 
described above, the minimum possible decay length L [estimated by substituting = y/s/2 
and y/F ~ A in eq. (0)] is typically of order 10 to 100 microns, far smaller than the resolution 
of the detectors. 

Neutralinos can be pair-produced in e + e~ collisions by Z exchange in the s-channel or 
by selectron exchange in the t-channel. In GMSB models, N± is always predominantly a bino 
and the ratio m^/m^ cannot be larger than (1.6, 1.25, 1.1) for n = (1,2,3). Therefore the 
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dominant contribution to N\N\ production always comes from e R exchange because of the 
relatively large e-e K -bino coupling. Indeed, the non-minimal GMSB models with n = 2, 3 tend 
to have a larger cross section for e + e~ — * AiAi because of relatively lighter mg R for fixed mft . 
In addition, in some models £r£r production (or just fifi) production is allowed at LEP2. 
However, we find that in the neutralino NLSP scenario, the individual slepton production cross 
sections are always smaller by at least a factor of 2 (and often much more), so we will concentrate 
first on N\N\ production as the discovery process. We also find that chargino pair production 
is never allowed at LEP2, and N1N2 is sometimes allowed but is always highly kinematically 
suppressed and therefore insignificant. This is easily understood due to the assumed bound 
N\ > 70 GeV and the rough relations rrift ~ ~ 2m ft , which hold since fi is relatively 
large in our models. 

The discovery process e + e~ — > NiN\ with each N± — > 7G, leads to events with two 
energetic photons and large missing energy.^ This signal has already been studied for LEP2 



in an earlier paper [30], but we will be able to extend these results. Also, ref. [3C] made no 
assumptions about model parameters; in taking into account the constraints inherent in the 
GMSB models we will be able to make some more concrete (and optimistic!) statements. The 
most important point to be made in this regard is that in the GMSB neutralino NLSP models 
of the class described in section II, the N\Ni production cross section can be bounded from 
below for a given m ft . This result is due to the facts that Ni is always gaugino-like and rrie R 
is bounded from above. The dominant t-channel exchange of e R therefore always ensures a 
substantial cross section. To illustrate this, we show in Figure [| the minimum and maximum 
cross sections for N± pair production obtained in these models, for ^fs = 161, 172, and 190 
GeV. This graph was prepared by an exhaustive scan of the model parameter space, varying 
A, x, tan/3, and sign(/i). The minimum cross sections are obtained for models with n = 1 and 
x not too large in which mg R /mft^ saturates its upper bound of about 1.6. 

The energy of each Ai is equal to the beam energy \fs/2, so that the photon energies in 
each event have a flat distribution, with 

-E'min < -B71 j -^72 ^ -E'max (32) 

where 

1 



-Emax.min = ^(\/s ± - 4m^) . (33) 

Therefore one can always impose a cut on soft photons, depending on an assumed lower bound 
on the mass of the Ni being searched for. In this paper, we will (motivated by |30|) take 



3 We prefer not to use the words "acoplanar photons" to refer to this signal, since the acoplanarity seems not 
to be a particularly useful discriminant against background. 
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Figure 2: The maximum and minimum cross sections for e + e~ — > N\N\ at ^fs = 161 GeV 
(short-dashed), 172 GeV (long-dashed), and 190 GeV (solid), as a function of the mass of 
the lightest neutralino. These bounds hold in the neutralino NLSP models within the GMSB 
framework described in section II. 



m 



> 70 GeV as a given, so that a cut 



1 



Ej > ~r(Vs - J s - (140 GeV) 2 ) 



(34) 



on soft photons can be applied without affecting the signal at all. The missing energy in each 
event is also bounded according to 2E m \ a < f < 2E max . 

The most important physics backgrounds for the 77^ signal are due to e + e~ — ► 771^17 
from diagrams with s-channel Z exchange (i = e,fj,,r) and t-channel W exchange (i = e only). 



These backgrounds were discussed in some detail in [3C], where it was shown that they can be 
efficiently eliminated using a cut on the invariant missing mass M : 2 NV = (p e + +p e - — p 7l — p 72 ) 2 - 
The Minv distribution of the signal tends to be broadly distributed and has most of its support 
for lower values than the background, which is strongly peaked at Mz but with a significant 
tail due to the t-channel contributions. In Figs. ||| and |] we show a comparison of the signal and 
background distributions for y^s = 172 and 190 GeV, respectively. In preparing these figures, 



we have applied detectability cuts 45] 



I cos&yl < 0.95 , 
(pr) 7 > 0.0325 ^/s, 

for each photon. The background distributions before and after the cut on soft photons (|34"|) were 



(35) 
(36) 
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Figure 3: Missing invariant mass distributions at y/s = 172 GeV for 77^ events which pass 
the cuts (35) and fl36|). The background from e + e~ —* 77^Fj is shown as the open (shaded) 



histogram before (after) the photon energy cut > 18 GeV from (|34|). The upper (lower) 
dashed line is the signal distribution from N±Ni production for representative models with 
= 75 (82) GeV as described in the text. 

computed at tree level using CompHEP p6[|, and are shown as the open and shaded histograms. 
Also shown are the distributions for the signals derived at tree-level from some sample models 
with m,fj i = 75, 82, and 90 GeV. These models were chosen with n = 1, tan /3 = 3, x = 0.1, and 
fi < 0. The shape (but not the magnitude) of the signal distributions is largely independent 
of these choices, for a fixed \fs and fn^ . This can be easily understood since the dominant 
kinematic features of these events are due to the isotropic decays N% — > 7G. Besides the cuts 
eqs. (^4|)-(^) we therefore propose to implement a cut on the missing invariant mass of 

5 GeV < Minv < 80 GeV. (37) 

The lower limit in ([}?]) is designed to remove detector backgrounds such as e + e~ — ► 77(7) 
with the third photon unobserved. We expect that the effect of initial state radiation will 
be to slightly reduce the total magnitude of both the signal and the background, while not 
qualitatively affecting the shapes of Minv distributions for our purposes. While significant 
quantitative changes are possible for Minv ^ m z, the cut eq. ( |37D makes these effects irrelevant 
here. 

At \fs = 172 GeV the total N\N\ production cross sections for the examples in Fig. || 
with =75 and 82 GeV are 195 fb and 41 fb respectively, before any cuts. After the cuts 
eqs. (j34|)-(|37|), the remaining background is less than 1 fb, while the signal for these models 
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Figure 4: Missing invariant mass distributions at y/s = 190 GeV for 77^ events which pass 
the cuts (35) and p6|). The background from e + e~ —* ^v^pi is shown as the open (shaded) 



histogram before (after) the photon energy cut > 16 GeV from (34). The dashed lines are 
the signal distributions from N±Ni production for representative models with = 75, 82, 
and 90 GeV (from top to bottom) as described in the text. 

is 128 fb and 30 fb. With 10 pb _1 per detector, this amounts to an expectation of roughly 5 
events and 1.2 events (summed over all four detectors). 

At y/s = 190 GeV the total cross sections in the models of Fig. § are (350, 180, 41) fb 
for = (75, 82, 90) GeV, respectively. After applying all of the cuts eqs. (j34"l)-(j3^), we find 
remaining signals of (193, 106, 25) fb. With 300 pb _1 per detector, this corresponds to about 
(58, 32, 7.5) events for each of the four detectors. The remaining background after these cuts 
is less than 1 fb. 

As shown by these examples, the efficiency for detecting signal events after cuts is quite 
high, while the physics background is essentially completely eliminated. The effects of the 
cuts on the backgrounds are shown in Table 1. In Fig. ||, we show the efficiency (defined as 
the fraction of signal events which pass all of our cuts divided by the total number of signal 
events) as a function of for various beam energies y/s = 161, 172, 185, 190, 195, and 200 
GeV. The efficiency decreases slightly with increasing beam energy, but always exceeds 50%. 
Also, we note that the efficiency increases slightly closer to threshold; this is because nearer 
threshold the Mjnv distribution for the signal becomes somewhat more sharply peaked with 
a smaller overlap with the cut region Mjnv > 80 GeV. The plotted efficiencies were found for 
a specific class of models with n = 1, tan/3 = 3, x = 0.1, /U < and varying A, but these 
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Figure 5: The fraction of 77^ events from neutralino pair production which pass the detector 
cuts ([35]) , ( |36| ) , and missing invariant mass cut ([37]) . The lines from top to bottom are for y/s = 
161, 172, 185, 190, 195 and 200 GeV. Here we have chosen a representative class of models as 
described in the text, but for fixed rrijj and y/s the efficiency is quite insensitive to variations 
in model parameters. 



results are very nearly model- independent for a fixed m pj , since the efficiencies depend mostly 
on the kinematics of the isotropic JVi decays. Therefore one may estimate the usable cross 
section (after all cuts) for future LEP2 runs by simply multiplying the total cross section by 
0.5. While the LEP runs at y/s = 161, 172 were limited by only having ~ 10 pb" 1 collected 
per experiment, the reach at y/s = 172 GeV extends even beyond =77 GeV in some 
parts of the GMSB model parameter space. However, considering the minimum cross section 
obtained in some models (see Fig. |2|), the exclusion capability for these runs is about > 72 
GeV, using a criterion of 5 total events (summed over all four detectors) after all cuts. This 
limit is quite comparable to what should be attainable with the present Tevatron data (and 
provides a posteriori justification for our assumption > 70 GeV). With 300 pb^ 1 and a 
discovery requirement of 5 events after cuts, the discovery reach should extend up to about 
m N 1 = — 4 GeV in future runs, based on an efficiency of 50%. If a discovery is made, the 



events with the largest photon energies can be used to find using eqs. (32) and (|33D . (The 
lower endpoint of the energy range will be contaminated with "fyvV background events.) 

We note that observation of a few 77^ events with Minv > 80 GeV could not be interpreted 
as unambiguous evidence for A7A7 production, since the background is comparable to or larger 
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than the signal, especially if our cut on soft photons ([34]) is not applied. Conversely any events 
with Minv < 80 GeV would be of great interest for the GMSB scenario, since the physics 
backgrounds for such events are quite negligible.^] 

When kinematically allowed, slepton pair production can add to the main NiNi discovery 
signal. Each of the production cross-sections for e R e R , p- R p, R , and fifi can even exceed 0.1 pb at 
y/s = 190 GeV, adding events l + £~^^Ifi to the 77^ signal we have just discussed. The leptons 
in these events should be softer than the photons, because the slepton production cross section 
cannot be significant unless mj—m^ is small. For some models with large tan (3, only f\f\ pair 
production with signal t + t~77^ can occur in addition to 77$. While the individual slepton 
pair production cross sections never exceed half of the N\N\ production cross section, the 
£ + ^~77^5 signal(s) should not have any significant backgrounds (especially considering that the 
photons are always energetic) . Therefore it should be kept in mind that slepton pair production 
can be an important component of the discovery signal even in the neutralino NLSP scenario. 
In the minimal model with n = 1, this can only occur if x is significantly greater than 0, so 
that observation of both 77^ and £ + £~^^]fi at LEP2 would exclude models with n = 1 and 
small relative mass splittings for the messenger fields. 

B. The slepton co-NLSP scenario 

In this subsection we consider the case that f%, e R , and fi R are lighter than all of the 
other superpartners, and are nearly degenerate in mass. As long as the conditions eq. ([ll]) are 
satisfied, then f\, e R , £l r act effectively as co-NLSPs, each decaying directly to the corresponding 
lepton plus gravitino. This situation can arise if n = 2, 3, or 4, and tan (3 < 8. For larger 
values of tan {3, the tau Yukawa coupling causes mixing between the left- and right-handed staus 
which is always sufficient to render lower than mp, R by more than the r mass. Restricting 
our attention to models with slepton NLSP masses between 50 and 100 GeV, we find that the 
allowed ranges for the parameter A are 

n = 2; 15 TeV < A < 42 TeV, or (38) 
n = 3; 11 TeV < A < 35 TeV, or (39) 
n = 4; 10 TeV < A < 28 TeV (40) 

in order for eq. to be satisfied. As mentioned in the Introduction, however, there are some 
models for which ra~ — — m T — rri£ can be up to a few GeV but In — > IG can still dominate 

4 We also note that in the higgsino LSP interpretation |4?| of the CDF event Jis]] , one could conceivably 
have 77-^ events from N2N2 production with the one-loop decay N2 — * 7-ZVi, but these would yield softer photons, 
and like the background would tend to have larger Minv than a GMSB signal. 
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if \/F is not too large and £r —* £N\ is not kinematically open. Those models will also act 
as slepton co-NLSP models. In any case, using eqs. (||) and (^5|), we find that the minimum 
possible decay length at LEP2 in the slepton co-NLSP scenario is about 10 microns for n = 2, 
and somewhat smaller for n = 3,4. Of course, values of x\ smaller than 1 will increase the 
decay length proportionally to l/(x\) 2 , and F$ < F would have the same effect. This means 
that the I — ► IG decay lengths can easily exceed minimum detector resolutions, providing a 
background-independent signal (2^, [34|] if the tracks of stable sleptons and/or their macroscopic 
decay lengths are observed by the LEP2 detectors. This would be spectacular confirmation of 
the GMSB scenario. 

The cross section for fi R jl R production at LEP2 (as a function of to^) is model-independent, 
since this process is mediated only by s-channel exchange of 7, Z. The e R e R pair production 
cross section has a contribution from i-channel neutralino exchange, but there is significant 
destructive interference with the s-channel 7, Z exchange graphs, especially near threshold. 
This means that a{e R e R ) is often much lower than a{jl R jl R ) in slepton co-NLSP GMSB mod- 
els. This is an important qualitative difference from the situation in neutralino LSP models as 
studied in , where exchange of a lighter neutralino in the t-channel typically ensures that 

a(e R e R ) > a(fl R jl R ). The destructive interference effect in slepton co-NLSP models is greater 
for larger values of n, corresponding to heavier Ni. We find that in our slepton co-NLSP model 
parameter space, a(e R e R ) < <r(p, R p, R ) always holds for raj more than about (10, 16, 20) GeV 
below the kinematic threshold of y/s/2, for models with n = (2, 3, 4). 

These features are illustrated in Figs. || and [7], which show production cross sections in 
e + e~ collisions at = 172 and 190 GeV for e R e R and jl R fl R . The results shown for e R e R are 
for a typical family of GMSB models with x = 0.1, tan/3 = 1.5, /x < and A varying, as a 
function of m^ R rn,e R - The three dashed lines in Fig. || and Fig. [7| are the e R e R cross sections 
in this class of models for n = 2, 3, 4, from top to bottom. (These models have lower bounds on 
m~ as indicated; these follow indirectly from a lower bound on the mass of the lightest Higgs 
boson, which we take to be m h < s'm 2 (/3 - a) 64 GeV @.) For a given n, the e R e R production 
cross sections in other models can be up to a factor of two smaller, but not much larger, than 
shown in Figs. || and [?]. The production cross section for f\f\ is always nearly equal to that for 
fl R £i R , owing to the small mixing required by low tan/3 and near-degeneracy of jl R and f\. For 
a given model, one finds that g{t\Ti) exceeds a(jj, R ji R ) by a few per cent. 

With the 10 pb _1 collected at \fs = 172 GeV it is apparent from Figure ^ that, modulo 
detector-dependent considerations regarding the rate of energy loss by ionization and tracking 
chamber capabilities, one should be able to put a useful exclusion on long-lived sleptons. Indeed, 
the DELPHI collaboration has recently analyzed data from runs with yfs < 161, 172 GeV, and 
found that long-lived right-handed smuons and staus with mass less than 65 GeV are excluded 
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Figure 6: Slepton pair production cross sections as a function of slepton mass, for e + e~ collisions 
at = 172 GeV. The solid curve is the (model-independent) production cross section for jl R ji R . 
The three dashed lines are the e R e R production cross sections for a family of models with varying 
A and fixed x = 0.1, fi < 0, tan/3 = 1.5, and, from top to bottom, n = 2,3,4. (The model 
dependence is due to the t-channel exchange of neutralinos.) The f\f\ production cross section 
in slepton co-NLSP models is always nearly equal to that of fi R fi R . 



at the 95% confidence level [49|. A somewhat more restrictive lower bound could presumably 
be obtained by combining the results of all four detectors. Likewise, searches for sleptons with 
decay lengths exceeding several centimeters could probably establish similar limits. Comparing 
with Figure |?], we see that in future runs with 300 pb _1 or more, it should be possible to 
exclude or discover long-lived sleptons with masses up to a few GeV of the kinematic limit. 
(Comparable constraints on long-lived sleptons from the present Tevatron data will probably 
be difficult to obtain 



In the remainder of this subsection, we concentrate on the more difficult situation that 
the finite decay lengths for the sleptons are too short to measure. The signals for slepton pair 
production are then e + e~$, fj, + n~$, or t + t~$ for e R e Rl jl R fl R , f\f\ respectively ^6|. As 
illustrated by the examples of Figures ^ and ^, the worst-case situation for a given mj will 
have a(e R e R ) <C a(fj, R jl R ), so we concentrate on JX R jl R as the discovery process. (In general, 
e R e R can be the dominant discovery process in slepton co-NLSP models only when both cross 
sections are large and discovery is relatively easy anyway.) In Figure ||, we show the model- 
independent production cross section for jl R fi R as a function of mp, R , for various beam energies. 
The past LEP runs at </s = 130-136, 161 and 172 GeV collected about 5.7 pb -1 , 10 pb" 1 and 
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Figure 7: As in Figure [|, but with i/i = 190 GeV. 

10 pb _1 per detector respectively. Slepton masses less than about 55 GeV would have resulted 
in several d. + t~f events for each detector at LEP130-136, with a detection efficiency probably 
well in excess of 50%, since the decays i — > IG will always result in energetic leptons. (See for 
example the analogous situation analyzed in [^] in the case of sleptons decaying to a lepton 
and light neutralino.) For the LEP161 and LEP172 runs, several events per detector could 
be expected for slepton masses up to perhaps 70 GeV, but with a background from W + W~ 
production with leptonic W decays, as discussed below. A precise determination of slepton 
mass exclusions from present data in the slepton co-NLSP case would involve detector-specific 
issues and will not be attempted here. 

Future runs with higher beam energy and much more data should be able to decisively 
probe a significant range of slepton masses. The lepton energies in each event have a flat 
distribution (before any cuts) with endpoints given by 

-Emiri < E e +,E e - < -Emax (41) 

where 

£max,min = j(Vs ± yjl - 4m|) . (42) 

Therefore the leptons from the signal events are quite energetic (especially in the critical case 
that rriji is near the kinematic threshold so that the production cross section is low), allowing 
one to choose a rather strong cut on the minimum lepton energy. To reduce an important 
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Figure 8: Smuon pair production cross sections as a function of mn R , for e + e collisions at 
(from left to right) y/a = 130, 136, 161, 172 (solid lines) and 185, 190, 195, and 200 GeV (dashed 
lines). 

component of the background as discussed below, it is necessary to impose such a lower bound 



on lepton energy (in contrast to the situation for slepton signals with a neutralino LSP [43, 44], 
where it is instead useful to impose an upper bound cut on lepton energies). We will somewhat 
arbitrarily take this to be 

Ei > 20 GeV, (43) 

which does not impact the signal at all for mj> 77.5 GeV for y/s = 190 GeV, but this can and 
should be adjusted depending on the signal and on the collider parameters. In the following 
discussion we also impose a detectability cut 

W\ < 2.5 (44) 

on the pseudorapidity of each lepton. 

The fraction of e + e - ^, // + /U _ ^ and t + t~^ signal events which pass these cuts is always 
quite high. Unfortunately, there are significant backgrounds to these signals which must be 
considered. First, one has e + e~ — * t + t~ with leptonic r decays which gives a background 
to l + £~]fi with £ = e,/U. The resulting leptons are always nearly back-to-back, so that an 
acoplanarity cut 

cos(j)(£ + £-) > -0.9 (45) 
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can essentially eliminate this background, while leaving the signal largely intact. (The angle 
4> is defined between the two leptons in the plane transverse to the beam axis.) Backgrounds 
from £ + £~j and £ + £~e + e~ production with the photon or e + e~ lost down the beampipe can 
be efficiently eliminated with a cut on the total missing transverse momentum |j43"|: 



j> T > 0.05^. (46) 

Significant backgrounds also arise from e + e~ — > Z(—> vV)^f*{-^ £ + £~) and e + e~ — > Z(—* 
uV)Z{-^ £ + £~). To eliminate them we impose a cut 

|Mi NV - M z \ > 10 GeV (47) 

on the missing invariant mass Minv — {p e + + P e - ~ Pp+ ~ Pp-) 2 m each event. Other minor 

i i f f 

backgrounds are present from higher order processes such as e + e~ — ► ez/W(— > ev). Their sizes 
before and after our cuts are given in Table 2, where we show that they are reduced to a 
negligible level. It should be noted that these cuts will also reduce any potentially dangerous 
interferences of these backgrounds with the main background we are about to discuss. 

When yfs > 2M\y, the largest physics background is due to VF-pair production followed 
by leptonic W decays: 

e + e~ -» W + W~ -> £+£-uu. (48) 

At y/s = 190 fb, this background amounts to 238 fb for each lepton flavor, before cuts. This is 
reduced to 140 fb after applying the cuts eqs. (03)-(H7]). The kinematics of the background (H8j) 
are similar to those of the slepton pair-production signal, especially if m^-is close to mw, and in 
particular the cut (|43| ) has only a very small effect. The situation is rather similar to the case of 



slepton pair production at LEP in the neutralino LSP scenario as studied in 44], but with 
£ — > £G taking the place of t — > £N\. Note that in the present situation, the near masslessness 
of the gravitino makes the decay £ — > £G even more kinematically similar to the Standard 
Model decay W — ► £v. In particular, the signal cannot be enhanced significantly by imposing 
an upper bound on the lepton energies, as it could be in the situation investigated in |43|, EJJ. 
However, we can still use the fact that the positively (negatively) charged leptons from W + W~ 
production are produced preferentially in the same direction as the positron (electron) beam. 
This polar angle asymmetry also unfortunately holds true for the e + e - ^ signal from e R e R 
production, although not as strongly as for the background. The ^ fj,~ ~Jf) and r + r~ f signals 
are fortunately symmetric with respect to 9 — > tt — 9 because e + e~ — > p, R fiii and e + e~ — > f\f\ 
production have only s-channel contributions. Therefore the signal/background ratio for fi R fi R 
and f%T\ (and to a lesser extent e R e R ) can be significantly enhanced by imposing a cut 

± cos 9 e ± > (49) 
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on the more energetic lepton in each event, as in [43, |44j]. (We always use the definition of 9 as 



the angle between the e beam momentum and the outgoing lepton momentum.) 

We must also consider a background contribution for e + e~$ and (but not for 

t + t~|5) from 

e + e" -> W + W~ -»■ ^r^i/F -► t + rvVW (50) 

Since r — » IvV has a branching fraction of about 0.18 for each of I = e, fj,, a rough estimate is 
that the additional £ + £~]ft backgrounds are each about 0.36 times the direct W + W~ — ► i^t^vv 
background. However, this is an overestimate, since the resulting leptons tend to be softer so 
that a large fraction of these events are eliminated by the minimum energy cut which 
was included for this reason. Before (after) the cuts eqs. fl43])-(f4"7|) this background contributes 
84 (18) fb. It should be noted that if some efficient tau tagging is possible, then some part 
of this background could be eliminated. However, this is a highly detector-dependent matter 
and we choose to simply consider the whole background. There is also a background from 
e + e~ — ► W + W~ — ► t + t~vV — > l + £~vVvv : but this is suppressed by the factor BR(t — ► ivV) 2 ~ 
0.03 and is greatly diminished further by the cut (f43|), and so can be safely neglected. In Table 2 
we summarize the dominant l + £~$ backgrounds showing the effects of the cuts described above. 



To evaluate the background from eq. ( |50|) we used our own Monte Carlo simulating the 
tau decay following each of the parent WW — * [itvv events generated using CompHEP and 
the BASES/SPRING package [52|. In doing this, we also took into account spin correlations 



and tau-decay anisotropies in the tau rest frame. We checked that neglecting the latter effects 
would have resulted in an 18% underestimate of this background after cuts (1-5) in Table 2 and 
in an only slightly flatter distribution in ± cos 0£± so that the underestimate is diminished to 
17% after cuts (1-6). The effect of taking into account spin correlations is not dramatic in this 
case because the pattern of distributions is dominated by a large boost of the tau in the lab 
frame. 

These considerations are illustrated in Figure [9|, which shows the distribution of ± cos 9g± 
for the pL R fiR and e R e R signals in a model with mj = 80 GeV. (The other model parameters, 
which do not affect the distribution from smuons, were chosen to be n = 3, A = 24 TeV, 
tan/3 = 1.5, x = 0.1, with /i = —460 GeV.) Also shown as the heavier solid line histogram is 
the distribution for background events from both ( [48|) and (50). The /i 4 ' ji~ 'f! component of 



the signal is much more promising than e + e~]ft, both because the total cross section is larger 
and because the polar angular distribution is less similar to the background. This is a quite 
general feature. In the case of the T + r~]fi signal, there is no background from (|50"|), so we 
also show in Figure |9| the distribution from ( |48| ) only, as the lighter solid line histogram. (The 
T + r~]fi signal is, however, subject to a significant and quite detector-dependent loss from r 
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Figure 9: The differential cross section da/d(± cos 9#±) for the more energetic lepton in each 
event at y/s = 190 GeV, after the cuts ([43|)-(ff7|). The long-dashed (short-dashed) line is the 
signal from e + e~ — > £ + £~ — > £ + £~fJ for muons (electrons), respectively. For this figure we 
have chosen a model with m~ = 80 GeV as described in the text. The lighter solid line 
histogram shows the background contribution from e + e~ — ► W + W~ — > £ + £~vv (not summed 
over lepton flavors), while the heavier solid line histogram includes also the background from 
e + e~ — > W + W~ — > £^t^vV — > £ + £~vVvv where £ = e or /i. 



identification efficiencies.) 

As these examples illustrate, we may concentrate on fj> R jl R production as the likely discovery 



process, at least in the pessimistic but common situation that a{e R e R ) < a(jl R jl R ). In Figure 10 
we compare the distributions for the polar angle of the more energetic muon from background 
and signal events, for various jl R masses. The W + W~ background should contain e + e~]fi, 
e + /i - ^ and e~/i + ^ events in equal amounts. Meanwhile the signal yields more fi + fj,~^J 
than e + e~ ]fi events after the cut eq. (pfij), and no e/x^ 1 events [except from f\T\ production with 
both taus decaying leptonically; this signal is again strongly diminished by the square of the 
leptonic branching fraction of the tau and by the energy cut eq. (|43|)]. The cross sections at 
yfa = 190 GeV for signal /i + >" $ events before cuts are (360, 238, 178, 121, 69, 25) fb for 
m VR = (60, 70, 75, 80, 85, 90) GeV, respectively. The corresponding amounts which pass all 
of the cuts including eq. ( |4G|) are (70, 60, 52, 39, 23, 8.5) fb respectively. This amounts to an 
expectation of (21, 18, 16, 12, 7, 2.5) signal ^ \i~ f, events per detector with 300 pb _1 . For 
comparison, the total background from eqs. (f48|), ( |50| ) and other sources which survives 

all of the cuts including eq. (fl9|) is 26 fb (see Table 2), or 7.8 events in 300 pb _1 . Combining 



the results of all four detector should therefore give a 5a discovery signal for mj^ up to 85 GeV. 
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Figure 10: The differential cross section da/d(± cos 0£±) for the more energetic muon in each 
event at = 190 GeV, after the cuts ([43|)-(|47|), The dashed lines are the signal from e + e~ — > 
UrP'r — > H + f° r ) from the top down, mp, R = 60, 70, 75, 80, 85, 90 GeV. The background 
distributions are shown as described in Figure ||. 

However, we conclude from Fig. [l(] that an unambiguous discovery may not be possible for 
models with m,f lR > 85 GeV, even with 300 pb _1 per experiment at -y/i = 190 GeV, because 
of the low rates and significant WW background^] 

In the delicate region where one is searching for smuons with m^ R > 80 GeV, it will help 
to increase the minimum muon energy cut eq. (^) and to impose an upper bound cut on muon 
energies, in view of eq. (fl2|), since this will reduce the W + W~ backgrounds somewhat. It 
may also be useful to adjust the polar angle cut (fl9|). The optimal cuts clearly depend in a 
non-trivial way on the beam energy, on the amount of integrated luminosity available, and on 
the masses of the sleptons being searched for, because of the low signal rates. If a discovery 



is made, the events with the highest muon energy eq. (43) should provide the best estimate of 
the smuon mass using (|42"D, if m^ R < 80 GeV. For mp iR J> 80 GeV, such a determination will 
be more difficult because of the lower rates and because the range of muon energies from the 
signal events is entirely covered by the background. 

It is important to note that when kinematically allowed, N%Ni will likely provide the 

clearest discovery signal. Despite > , mg R , rrifi R in these models, we find that a(NiN\) 

J In the above analysis we did not take into account initial state radiation effects, which diminish total cross 
sections for both the signal and the backgrounds by up to 10%, and produce slight changes in the shapes of 
distributions. We do not expect this to significantly affect our conclusions. 
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can be as much as 3 times larger than a(jl R p, R ) in the slepton co-NLSP scenario. Each Ni 
can decay as N\ —>■ £r£ — * £ + £~G or N% — * f\T — > t + t~G. The first lepton emitted in the 
decay chain will often be very soft (if — mj R or - — rn T is small), but this does 
not degrade the signal. The reason is that because of the Majorana nature of N\, the charges 
and flavors of the two most energetic leptons in the event (which carry most of the visible 
energy) are uncorrelated with each other. Therefore one has the rather unique signature of 
£ + £' + {£-£'-)$ and £~ £'~ {£ + £'+)$ and £+ £'~ {£~ £'+)$ events in the ratio 1:1:2, with the leptons 
in parentheses being much softer, perhaps even too soft to detect. Here £ and £' can each be e, 
//, or r, independently. Even if both soft leptons are not detected, the presence of two energetic 
same-charge leptons with large missing energy should be an unmistakable discovery signal for 
the first two of these signatures. (The analog of this signal in the special case of a stau NLSP 
scenario was recently discussed in p7| .) The N\Ni production cross section in slepton co-NLSP 
models is still bounded from below as in Figure || (as a function of rrifi ), but can also be much 
larger. We therefore find that the NiNi channel will generally provide the largest discovery 
signal in slepton co-NLSP models if is more than about 5 GeV below threshold. This is 
particularly likely to happen in models with n = 2, for which Ni tends to be not much heavier 
than the sleptons. It is also especially probable when cr(e R e R ) and cr(jl R fl R ) are also both large, 
so that identifying an excess of e + e~$ and events above the WW backgrounds will 

also not be difficult. In that case the slepton pair production and N\N\ pair production signals 
should provide strong confirmation of each other. 



C. The stau NLSP scenario 



In this section we consider the case that the lightest stau is the NLSP and all supersym- 
metric decay chains terminate in f\ — * tG. Restricting our attention to f\ masses between 50 
and 100 GeV, we find that the condition eq. ( |l0| ) can be satisfied in the GMSB model parameter 
space of section II only if 



n 
n 
n 
n 



35 TeV < A < 120 TeV 

18 TeV < A < 80 TeV 

12 TeV < A < 70 TeV 

10 TeV < A < 60 TeV 



tan (3 > 18, or 
tan j3 > 6, or 
tan j3 > 5, or 
tan (5 > 4. 



(51) 
(52) 
(53) 
(54) 



These requirements are necessary but not sufficient for the stau NLSP scenario; indeed, as we 
have already mentioned, £r — ► IG for i = e,fx can dominate even if eq. (10) holds, provided 
that y/F is not too large and £r — * £N\ is not open. 



As in the slepton co-NLSP scenario of the previous subsection, the discovery prospects are 



27 



clearest if the decay length f\ — * tG is macroscopic, so that the discovery signal consists of 
tracks from a heavy fx and/or kinks due to f\ decays which can be directly observed in the 
detector |34|] . In this case, there are no significant physics backgrounds, and the discovery 
potential is limited only by the total production cross section, the integrated luminosity, and 
the capabilities of the detectors. In the GMSB models with a stau NLSP, the cross section for 
e + e~ — ► f\f\ as a function of is given to a good approximation by the curves in Figure || 
with the horizontal axis now interpreted as m^, and in particular is not very model-dependent. 
The mixing in the stau (mass) 2 matrix does provide for a small reduction in the f\f\ cross 
section compared to that shown, but we have checked that this reduction is at most about 10% 
in our models. As can therefore be seen from Figure ||, future LEP2 runs should be able to 
discover a long-lived stau NLSP with mass up to close to the kinematic limit, given 300 pb _1 
or more. 

In the following, we therefore concentrate on the more difficult possibility that the f\ decay 
length is too small to be directly observed. All events will then have an energetic t + t~ pair 
and large missing energy. In the large tan/3 limit, one has <C mg R ,m^ R ,m 1 ^ i in these 
models, so that the only discovery signal at LEP2 is t + t~^. This must be compared to a 
background from W + W~ with W — > rv decays, as given in the previous section. Since the 
signal events feature E T with a flat distribution as in eqs. and (42), the tau decay length 



of ~ 90 fj,m xE T /m T (roughly 1 to 4 millimeters) may allow for fairly efficient tagging of non- 
leptonic r decays. Of course, the W + W~ background produces taus with a similar energy 
distribution. Just as discussed in the previous section, the f\f\ signal is symmetric with respect 



to 9 — > 7r — 6, while the background is not. Modulo the tau identification problem, Figure 10 
gives an indication of the polar angle distributions of signal and background for t + t~^ events. 
Extraction of a f\f\ signal from the background will be considerably more difficult than would 
be the case for jl R jl R signal in the slepton co-NLSP scenario, and an estimate of the reach will 
depend quite sensitively on detector capabilities. If the taus are not tagged, one possibility is 
to look for purely hadronic states with very large missing energy to avoid WW contamination. 

For lower tan/3, the production of e R e R , fL R fl R , and/or NxNx can also be kinematically 
allowed. For a given model, the p, R p, R cross section is always smaller than that for f\f\ pro- 
duction, due simply to the kinematic suppression associated with m^ fl > + m T . It can 
of course be read off of Figure |8L as before. The e R e R production cross section can be either 
smaller or larger, due to the interference between graphs with t-channel exchange of neutralinos 
and s-channel exchange of 7, Z. If tan/3 does not exceed about 30, these cross sections can 
add to the signal. Now because of the decays £r — ► It fx-, the signal from IrIr production 
is t + t~ (l + £~T + T~)fl!, where the leptons in parentheses are much softer. If the leptons are 
energetic enough to be identified, this signal should have very low backgrounds. It is also pos- 
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sible that both £r — > f\rl and £r — > £G have a significant branching fraction, if the former 
decay is sufficiently kinematically suppressed and £r — ► IN\ is not open. The IrIr produc- 
tion could then lead to the additional signatures £ + t + (£~t-)$, fr~ {£ + t + )$, £ + t~(£-t+)$, 
£~t + (£ + t~)$, and, as in the slepton co-NLSP scenario, £ + £~$. 

If N\N\ production is allowed it can provide the dominant signal. In this case, each Aq can 
decay predominantly to rf\ and then to t + t~G. The final taus from each f\ decay will combine 
to carry most of the visible energy in each event, and their charges are uncorrelated because of 
the Majorana nature of Aq. This provides for the striking signatures 
t~t~ (t + t + )$ and t + t~ (t + t~)$ in the ratio 1:1:2, with the parentheses denoting soft particles 
as before. As in the case of slepton co-NLSP models, we find that the AqAq cross section can 
be up to 3 times larger than that of t\t\, despite the requirement mrifj > + rn T . When 
tan /3 is not too large, this signal can be the most visible one at LEP2 for stau NLSP models. 
If the decays Aq — > ££r are also allowed for £ = e,fi, one may obtain the same signatures but 
with two or four additional soft leptons from the cascade decays of Aq through £r. In our 
stau NLSP models we find that the AqAq production cross section as a function of is still 
bounded from below (but can be up to 50% smaller than indicated in Figure ^ in some cases). 
If rrtjy — 5 GeV, there should be at least a few events with very energetic same-charge 

taus and a pair of softer taus with the opposite charge, if more than 300 pb _1 is obtained. This 
is again especially likely to be the discovery signal in models with smaller n and tan (3 not too 
large. 

D. The neutralino-stau co-NLSP scenario 

Finally we consider the case that the lightest neutralino and the lighter stau act effectively 
as co-NLSPs. This scenario will occur if the conditions of eq. ( |l^ ) are satisfied. While it might 
seem at first that requiring m ^ and to be nearly degenerate requires some fine-tuning, we 
find that the region of parameter space where this occurs is actually quite significant. Conditions 
which are necessary for a viable model (with ^ 70 GeV) in the neutralino-stau co-NLSP 
scenario are: 

n = 1; 40 TeV < A < 80 TeV; 15 < tan/3 < 40, or (55) 
n = 2; 20 TeV < A< 40 TeV; tan/3 < 20, or (56) 
n = 3; 17 TeV < A < 25 TeV; tan/3 < 12. (57) 

Because Aq and f\ masses are not very different in this scenario, strict exclusion or discovery 
should be rather straightforward for a given NLSP mass. This is easily understood in terms of 



37| t+t+{t-t- 



and 
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the other scenarios we have already studied. We find that in neutralino-stau co-NLSP models, 
the photino content of N% is bounded from below by «i 7 > 0.35. Comparing eqs. (Q) and (||) we 
therefore conclude that if the iVi decay always occurs outside of the detector, then the ft decay 
must also take place over a typically macroscopic distance. If N\N\ production is kinematically 
allowed, then as in section III. A one has a signal 77^ with about a 50% detection efficiency 
after cuts, as long as the N\ decay length is not too long. Conversely, if the decay lengths 
are long, then t\t\ production leads to a background-free signal with tracks of quasi-stable fi 
observed directly in the detector using the rate of energy loss and/or kinks due to slow decays 
fi — * tG. In an intermediate regime, one should see both types of events. Note that in any case 
one need not rely on identifying a r + r~|5 (or £ + £~f!) signal against the W + W~ background to 
effect discovery. In neutralino-stau co-NLSP models, we find that the N\N\ production cross 
sections at LEP2 are bounded from below as in Figure ^. Likewise Figure |8| can be used to 
estimate the f±fi production cross section as a function of m^. This overestimates the actual 
fifi rate by no more than 10% because of stau mixing effects. When and are more 
than about 5 GeV below threshold, it cannot be possible for both of these signals to elude 
detection in future LEP2 runs. 



IV. Concluding remarks 



In this paper we have studied the implications of gauge-mediated supersymmetry breaking 
models for LEP2 physics. There are four main scenarios for the effective NLSP(s), each with 
its own predictions for the possible discovery signals. These possibilities also depend on the 
unknown NLSP decay length, and are summarized in Table 3. 

In many cases, strict exclusion (or discovery!) for a given NLSP mass can be assured given 
a sufficient beam energy and integrated luminosity. For example, we found that a(N\Ni) is 
bounded from below (for a given ) in the neutralino NLSP scenario models discussed in 
Section II. We analyzed the backgrounds and found that they can be eliminated with cuts which 
retain at least 50% of the signal events. This should guarantee discovery of N\ with masses up 
to a few GeV of the kinematic limit in future LEP2 runs, provided that the N\ decay length is 
not too long. As a caveat, we must note that with variations of the model choices we have made, 
it is quite possible to find smaller cross sections for N\N\ production. For example, there could 
be additional corrections to the Higgs soft (mass) 2 parameters with an indirect result of smaller 
values for \^\ [ B4j . Similarly, we have investigated some models with unequal non-zero values 
for the messenger multiplicities J2i n i{^)i Y^i n 2(i) and Y^i n 3^) [instead of eq. (|20|)1, and found 
that the requirements of correct electroweak symmetry breaking can and do lead to smaller 
values of \fi\/M2 in viable neutralino NLSP models. In both cases, one finds that Ni can have 
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a large higgsino content, with therefore an arbitrarily small production cross section for a given 
and a very long decay length for N\ — > 7G. 

In the slepton co-NLSP scenario, we found that if decays are prompt, the signal 
from p, R fl R production is likely to be the discovery mode. In particular, we found that the 
e + e~ \fi signal from e R e R production can be highly suppressed by interference effects, and has a 
comparatively unfavorable polar angle distribution. We also found that it is necessary to employ 
a different cut on the lepton energies than would be used in the neutralino LSP scenario to 
separate the signal from the WW backgrounds, because of the kinematic characteristics of the 
decay £ — > IG with a nearly massless gravitino. The pair production of N\N\ can also lead 
to spectacular signatures involving two energetic leptons with the same charge (and two softer 
leptons with the opposite charge), in both the slepton co-NLSP and stau NLSP scenarios. Long 
slepton lifetimes should lead to a relatively easy discovery from observation of heavy charged 
particle tracks and/or decay kinks. One of the more difficult scenarios involves a stau which is 
much lighter than all of the other superpartners and which decays promptly into tG. In this 
case the only discovery signal is t + t~ 'f, with a significant background from WW production. 

It should be noted that the signals we have studied are considerably more general than the 
models outlined in section II. The same scenarios and qualitative features of the signals arise in 
a much larger class of GMSB models with, for example, different numbers of messenger fields 
and/or widely different messenger scales. In most cases the discovery of a GMSB signal will 
be readily distinguishable from the predictions of models with a neutralino LSP. Therefore the 
variety of different signal possibilities points to the exciting prospect of simultaneously discov- 
ering supersymmetry and uncovering some of the most prominent features of the mechanism of 
supersymmetry breaking. 
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Cross section (fb) for E 


771/j 


Vi after 




i=e,ii 


T 




(GeV) 


cuts: (1-2) 


(1-3) 


(1-4) 


172 


131 


65.6 


0.79 


190 


102 


70.0 


0.80 



Table 1: The ^vv background to the 77^ signal at LEP172 and LEP190, after the following 
cuts (described in the text): (1) |cos6> 7 | < 0.95, (2) (p T ) 7 > 0.0325^5, (3) £ 7 > \{yfs - 
y 1 s — (140 GeV) 2 ), (4) Mjnv < 80 GeV. The missing invariant mass cut (4) clearly reduces the 
main background to a negligible level. The additional cut Minv > 5 GeV is needed to eliminate 
the background from 77(7), where (7) is lost in either the detector or the beam pipe. 



Background {£ = e or //) 

a) W + W~ -> e + £-v e V e 

b) W+W- -► ^^(-^ e*v T vi)v t v T 

c) E 7*H^")^H^0 

i=e,/j,,T 

d) £ e + e _ Z(^ i/jFj) 

i=e,ix,r 

(other than ZZ, 7*Z contribs.) 

e) Y.±^eW^(^e^u e ) 

f) e 2h^ + r)^^) 



Cross section (fb) at LEP190 after 
cuts: (1-3) (1-4) (1-5) (1-6) 



163 


141 


140 


22.3 


57.2 


42.5 


17.7 


3.03 


31.2 


2.97 


0.83 


0.40 


14.8 


0.93 


0.56 


0.13 


13.9 


11.4 


7.69 


1.96 


4.92 


0.27 


0.27 


0.13 



Table 2: The dominant backgrounds to dilepton signals at y/s = 190 GeV after the following 
cuts (described in the text): (1) \rj e \ < 2.5, (2) i> T > 0.05^5, (3) cos(f>(£ + £~) > -0.9, (4) 
\M my -M z \ > 10 GeV, (5) E e > 20 GeV, (6) ±cos0 € ± > 0. Other channels [e.g. uVZ{-^ £ + £~); 
H ± i/ fl W^ : (—> [M^Vfj); contributions to fi + fi~Z from processes other than j*Z and ZZ production; 
multiperipheral diagrams; etc.] produce additional backgrounds at the level of 1 fb or less, before 
cuts. Cuts (1-2) are needed to reduce backgrounds from e + e~('y), e + e~ (£ + £~), where (x) means 
x is lost in the beam pipe, as well as to ensure final state detection. Cut (3) is necessary to 
eliminate the background from t + {-^ £^ ~vv)t~ "(— ► £~vv) and is also quite effective in reducing 
a), b), d) and e) [by about 24%, 23%, 18% and 15% of the corresponding respective amounts 
after cuts (1-2)]. Note that the processes d) and e) are backgrounds only for the e + e~ ~$ signal. 
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Sparticlc 

Scenario Production Signal Comments 


ncutralino NLSP < 


N\ — > 7G decays are 

displaced 7s 


j prompt 

j within the detector 


P+p- < 

Z R Z R 


ii£ + £r$ , , , 

(as above) 

+ displaced 7s 


slepton co-NLSP £ + Jr R 

1 


' « f 

£—> £G decay kinks £ — » £G decays are < 
charged £ tracks [ 


prompt 

within the detector 
outside the detector 


N.N, 


' £+£' + {£-£'-)$, 
£-£'~(£+£' + ) fi, and {££') leptons are soft or undetected 
£+£'-(£- £' + )f 


stau NLSP fi+ff < 


fi — > rG decay kinks fi — > tG decays are < 
charged fi tracks 


j prompt 

within the detector 
[ outside the detector 


mm 


' r + r +/ T - r -j|i (tt) leptons are soft or undetected; 

t~t~ (t + t + )$, and possibly with 2 or 4 additional soft 
v leptons (e or /xj it iVi — ► t-tjj is open 


£+£^ T+T-(£+£-T+T-)f (£+£-t+t~) leptons 


are soft or undetected 


neutralino-stau N\Ni , 
co-NLSP fi+ff 


77^, T+T-f 

displaced 7 s and NLSP decays are < 
n — > rG decay kmks 

charged n tracks v 


Drompt 

within the detector 
outside the detector 



Table 3: The possible signatures at LEP2 in the four different NLSP scenarios. The notation 

(1) "prompt", (2) "within the detector", and (3) "outside the detector" refers to a NLSP decay 
such that the decay vertex is (1) close to the interaction region and not measurably displaced, 

(2) possibly resolvable with a detector, and (3) well outside the detector. In the neutralino 
NLSP and slepton co-NLSP cases, £ stands for e, /x, or r; in the stau NLSP case £ stands for e 
or fi. (*) In the slepton co-NLSP case, the particular £ + £~ Ifi signature which is most likely to 
be observable is as explained in Section III.B. 
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